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ABSTRACT 


Analyses are performed on 170 wind profiles which were obtained 
at the Suffield Experimental Station during the 1959-1967 period. 
In addition, 34 wind profiles obtained at Kerang, Australia, are 
considered in the analyses. The data consisted of measurements of wind 
speed and temperature at a number of levels up to 96 meters. A com- 
parison of the results obtained through the use of the power law, 
logarithmic law and log-linear law, is made. It is found that the log- 
linear law represents the data in the near-adiabatic conditions more 
accurately than does either the logarithmic or the power law. However, 
under strong inversions.the power law gives a better representation. 

It is shown that the wind profile data used in the analyses 


are represented well by a modified log-linear law in the form, 
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This modified wind profile law is compared with the Monin-Obukhov 
log-linear law and Brooks log-power law. It is found that the differ- 
ences among these three laws are not statistically significant. 

A number of results obtained in this study are found to be in 
agreement with those found by an independent investigation at the 
Suffield Experimental Station. The results from this investigation 
are also compared with those obtained at the Brookhaven National 


Laboratory and at a number of other locations. 
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CHAPTER I 
OUTLINE OF THE PROBLEM 


1.5 Introduction 

The air motion in the free atmosphere is very nearly horizontal 
and predominantly determined by a quasi-balance between the horizontal 
pressure force and the horizontal Coriolis force. This quasi-horizontal 
air motion is modified considerably in the surface boundary layer which 
extends from the earth's surface to a height of approximately 100 meters. 
In this layer the wind is subjected to a considerable frictional force 
whose magnitude varies primarily with the nature of the surface and 
with the temperature distribution in the vertical. Moreover, the air 
motion in this layer is, normally, turbulent which causes a continuous 
diffusion of properties from one region to another. 

With the development of the theory of atmospheric turbulence 
at the beginning of the twentieth century it became necessary to 
consider the vertical structure of wind and temperature near the ground. 
Investigations showed that turbulence was related to a set of external 
parameters which had to be determined from the vertical wind profiles 
measured under various stability conditions. Furthermore, information 
on the variation of wind with height was required in many practical 
problems. Construction of high structures, design of airplanes and 
launching of satellites are only a few examples of where the wind 
profile information was utilized. With the rapid technological advance- 
ments in the past few years, accurate determinations of the wind 


profiles near the earth's surface became essential. In many applications 
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the variation of wind speed with height has to be expressed mathe- 
matically. No definitive law describing the variation of wind with 
height in the surface boundary layer exists as yet. However, several 
wind laws, describing the wind variation with height under various 
stability conditions, have been proposed. It is the aim of this 
thesis to investigate how well experimentally measured wind profiles 
can be expressed by the existing wind laws, particularly those which 


are most commonly used in theory and practice. 


1.2 The Problem 

Investigations by Frost (1948), Panofsky (1963), Hansen (1966) 
and others have shown that various theoretically and empirically 
derived equations, describing the variations of wind with height, 
fit the actual observations with various degrees of accuracy. The 
three main approaches taken to find a general wind profile equation 
were the logarithmic law approach, the power law approach and more 
recently the log-linear law approach based on the Monin-Obukhov 
similarity theory. Several sets of wind profile data have been ob- 
tained and analysed in the past. Most of the early data examined 
a relatively shallow layer of the atmosphere, generally no higher 
than a few meters. In most cases the early data were fitted to 
either the power or logarithmic law. More recently, with the construc- 
tion of meteorological towers, wind profile measurements have been ex- 
tended to a height of 100 meters and in a few instances even higher. 


During the last few years profile laws based’ on the similarity theory 


«S- 
-sitem bseesrqxs od of etd Jogtod diiw besge baiw to notisixsy of 
diitw bitw to nolssizey edd guidiroesb wal avisiniteb off .ylisorsem 
[stevee ,vevewoH .i9y 26) eteixe rays yrabnuod s7stre sit ni tdgten 


auolttev tabsu tdgisd d3itw notselisy batw dj gniditoasb ,ewsl baiw 





aid lo mis edt ef 41 .bs2qotq.meed veri .enotsibnos yiilidsjs 


aslitozq briw baaessm yiis jnatuitegqxs. I lawewod eisgitesvat o3 elesi2 


| 
| 


dotdw seodt ylrsIuottxsq ,2wsl briw gatitaixe off yd beeestqxs ad mB ~ 


.sotjon1g bas yioon ot bsev yinemmos Jom ete . 

meldoxd sit S$. “J 

(2001) meensh , (E901) wletonsd ,(BACI) JaosT yd anotssghieeval ; 
ylisottigqns bas ylisotss10sd3 evotisy jsit awode sved atedse bre / 
-iigted daiw bniw Io enoltisitsv sd3 gnidizossb .emotseups bevitzeb 

esiT ,.¥yos81u995 to vaatnwah evotrsev djiw snotzsvisedo Isutos eid 322 
noltisupa s{iiorq batw [sisneg s bait o3 nsasi asdosorggs olan som? 
stom bos dosorgqs wsl rswog ont dosetqgs wsl olmdtixsgol sf2 stew 

. 


vordstudO-ninoM edt no beasd dosorqgs wel zssaiti-gol sd32 ¢laneos2 - 
-do asad svsed aasb s{tioxrg bniw to etsa Isisvs2 .ytosH? yi frelimte 
beatmaxs stab yixs9 siz io JeoM .teaq sd3 at boayisas bas beniss 
weriginl on yifszensg .sreiqeomts ef io asysi wolisde yleviaeler ® 


03 bs33ii sxsw sisb yliss od2 e9282 Jeom al 819390 wel 8 i 
¥y 


carseat date sacha ed0M “wel eee! to a 


4 how 

. + voy Ue - od ; - : cae _ i 

8 19 Twessm erie 1g boiw .eiswe iyi an Rees — 
: : ae Pe ir tars [ bo Do ‘ os 

eet ; : - ‘ : 
. atts ; ~ pace i re is 
¥ 4 A " ba ne . ? : 7 
2” «ATS r ” yi ' Vth <s 


v ve - 


: 


= FN 


fn ae “a 





=o. 
have gained increasingly more acceptance. Panofsky and Prasad (1965) 


write: 


"The Monin-Obukhov similarity theory has made 
it possible to treat interaction of mechanical 
turbulence and heat convection in a meaningful 
manner. Thus, it is now generally accepted in the 
description of wind profiles near the ground". 


In the analyses presented in this study the wind profile data 
extending from surface to heights up to 96 meters were used to in- 
vestigate which of the wind profile laws represented the actual ob- 
servations most accurately, and under what stability conditions. 
Calculations were also made to determine the parameters which appear 
in different wind profile laws. A statistical approach was used to 
establish which law was best for a particular situation. The sample 
of more than 200 wind profiles was.considered large enough to warrant 
such an approach. It was hoped that this investigation of the vari- 
ation of wind speed with height would add to the understanding of the 


atmospheric processes in the lowest layers of the atmosphere. 


1.3. Instrumentation and Data 

Various diffusion experiments conducted by the Department of 
National Defense at the Suffield Experimental Station, Alberta, over 
the last two decades required precise measurements of the wind 
profiles. The chief source of the wind profile data was the 96-m 
micrometeorological tower located at the station. Micrometeor- 


ological sensors were located at different heights on the tower. 
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The wind speed measurements were made with Sheppard-Casella 3-cup 
contact anemometers connected electrically to recording counters which 
were photographed at one-minute intervals. The sampling periods 
varied depending on the experimental requirements; however, the most 
frequent sampling period was ten minutes. The anemometers were placed 
at heights of 2, 8, 16, 32, 48, 64, 80, 92 and occasionally at 96m 
above ground level. The temperatures were measured at heights of 
0.5, 1, 4, 92 and occasionally at 13 m. For a number of wind profiles, 
instead of the temperature measurement at one-meter above ground, 
the temperature measurement at two-meters above ground was obtained. 
Thus, the temperature gradients for layers 4-0.5, 13-1 or 13-2 and 
92-1 or 92-2 m were obtained. 

Besides the 96-m tower two movable 16-m towers, with anemometers 
placed at 0.5, 1, 2, 4, 8 and 16 m, were used to measure the variation 
of wind with height during some of the diffusion experiments. Temper- 
atures on the 16-m towers were measured at 0.5, 4 and occasionally at 
8 m above ground. 

The site of the 96-m tower was gently rolling prairie grass 
with no large obstructions within two miles. The 16-m towers were, 
generally, placed a few hundred feet from the 96-m tower. During the 
summer months the surface in the vicinity of the towers consisted of 
Sparse grass six to ten centimeters high with very sparse stalks up 
to 30 cm. During the winter months the shorter grass was covered with 
snow and very sparse stalks up to 20 cm high protruded through the 


moderately smooth snow surface. 
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During the 1959-1967 period 170 wind profiles, together with 
the corresponding temperature gradients, were obtained. The wind 
profile measurements were taken at different times during the year 
and under various stability conditions. A summary of the main 
features of the Suffield data is given in Table 1.1. Occasional 
malfunctions of the instruments caused some omissions in the data; 
however, the number of the missing observations was small in compar- 
ison to the total number of observations. 

The temperature gradients (AT), which were represented as 
the temperature differences, in degrees Fahrenheit, between the top 
and the bottom of the various layers, varied from -7.0 to 3.8 in the 
0.5-4%m layer; from -1.0 to 7.2°in the 1-13 m layer; and from -4.6 
to 24.5 in the 1-92 m layer. In order that the stability conditions 
in the three individual layers could be compared, the temperature 
difference corresponding to the adiabatic lapse rate was calculated 
for each layer. It was found that the temperature difference of -0.1, 
-0.2 and -1.6 F would produce approximately adiabatic lapse rate in 
the. 3,5, 12 and 91 m thick layers respectively. The following 
stability conditions were then defined: inversion, isothermal, stable, 
adiabatic (neutral) and super-adiabatic. The limits for these 
conditions together with the number of observations in each interval 
are given in Table 1.2. 

In addition to the Suffield data a sample of 34 wind profiles 
obtained at Kerang, Australia by Swinbank (1964) was analysed. The 


Kerang data were selected primarily as a control sample because, 
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Table 1.1 


Summary of the main features of the Suffield Data 





Profile Number 

Characteristic Interval of Obs. 
Jan - Feb - Mar 50 
Time of Apr - May - Jun 2) 
the Year (Month) Jul - Aug - Sep 9 
Oct - Nov - Dec 106 
0900 - 1300 34 
Time of 1301 - 1700 60 
Day (MST) 1701 - 2100 66 
2101 - 2400 10 
001 - 090 4 
Wind Direction 091 - 180 19 
(Deg. True) 1315-2270 122 
271. -~ 360 20 
Missing =) 
Hh 2,00 116 
Height of SOs ns 50 8 
Profile (m) 50 > her 16 5 
hae weO 41 
eto KG) Ly 
Sampling t = 10 LOL 
Period (Min) B30 aster LO ale) 
tie 50 37 
Th 32 50 
Surface 5021 e102 35 
Temperature (F) 702 T> 50 61 
aU) 16 
Missing 8 
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according to the literature, it is one of the most accurate sets 

of observations available at the present time. The heights of the 
anemometers at Kerang were identical to those used in the 0.5 - 16m 
layer at Suffield. Thus a direct comparison of the results between 
the two locations was possible. Unfortunately, the comparison was 
possible only for the ee cases because the Kerang data were 


limited to unstable conditions. 


Table 1.2 


Definitions of stability conditions and number of observations in each 
interval (Suffield data) 


Stability orl) 3 - To or T99 = T, 
condition No. No. No. 
Limits of Limits of Limits of 
(F) Obs.* (F) Obs. (F) Obs .** 
Inversion 70.0 19 >0.0 36 iO m2. 64 
Isothermal 0.0 2 0.0 8 = (2. 7 
COLO aL to 0.2 
Stable - - -0O.1 2 -1.3 7 
to -0.3 
Neutral =O) es 3 -0.4 3 -1.8 6 
to -0.1 Coe Or to -1.4 
Super- 
Va BB Ney: Nom ian eo gal O Ie Ly, <-0.4 14 <-1.8 cel 
See ee ee ee a a 
* Profiles extending to a height h < 16m 
** Profiles extending to a height h 7 80m 
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CHAPTER II 
TURBULENCE AND STABILITY PARAMETERS 


2.1 The Reynolds Number 

Until the beginning of the nineteenth century little was 
known about the laws of internal friction of fluids. Systematic in- 
vestigations of the flow phenomena in pipes and channels began with 
Osborne Reynolds in the middle of the nineteenth century. Because the 
fluid flow in straight pipes and channels resembles in many ways the 
motion of air near the surface of the earth a large part of the 
present wind profile theory for the surface boundary layer was develop- 
ed from these investigations. 

Reynolds (1883) conducted a series of experiments on flow in 
long straight pipes. By introducing colouring matter into water 
flowing through glass pipes he demonstrated the existence of two 
widely different modes of motion. The motion was laminar in character 
at low velocities and turbulent (or eddying) at high velocities. 

The change from laminar to turbulent motion took place suddenly at 

a critical velocity which was directly proportional to the kinematic 
viscosity (vy) of the fluid and inversely proportional to the 

diameter (d) of the pipe. Reynolds concluded that the transition from 
laminar to turbulent flow depends only on the value of a dimension- 


less ratio given by: 
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Here U is the average velocity and v = ufo , where uy is 
the dynamic viscosity and 9 is the density. This ratio became 
known as the Reynolds number and is a fundamental quantity in fluid 
mechanics. 

The Reynolds number has been the subject of many subsequent 
investigations both practical and theoretical. It was found that 
in a viscous flow the viscous forces can be represented as uu/L? 
and the inertia forces can be characterized as ore By where L 
is any characteristic length. The ratio of the inertia to the viscous 
forces is identified as the Reynolds number, giving the relative 


importance of the two forces. 


202 The Richardson Number 

While the Reynolds number was particularly useful in defining 
the transition from laminar to turbulent flow in pipes, where a density 
gradient is absent, it proved to be less useful in the atmosphere where 
a density gradient is generally present. It was found that the exist- 
ence of the buoyancy forces near the earth's surface has a major effect 
on the atmospheric turbulence. According to Businger (1955), the 
important distinction in the surface boundary layer is not between the 
laminar and the turbulent flow but rather between turbulence caused 
by mechanical friction and turbulence caused by convection. For 
example, in the neutral atmosphere turbulence near the earth's surface 
is entirely due to friction. As soon as a heat flux exists turbulence 


will either increase or decrease through the liberation or absorption 
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of convective energy. 

Richardson (1920) postulated that the ratio of the rate at 
which buoyancy forces extract energy from turbulence to the rate at 
which it is supplied by the shear stress determines whether turbulence 
will increase or decrease, This ratio became known as the flux 
Richardson number (Re) and is a fundamental quantity in micrometeor- 


ology. 


Lumley and Panofsky (1964) write the flux Richardson number 


as: 


R, e mY (2.2) 
Opeeuwhe( 5U/A9z) 


where the quantities wO and uw can be defined in terms of the 


exchange coefficients for heat (Ky) and for momentum (Ky) as: 


~ Zn) 
_- a= We ands ire = bl ( 


where 6 is the mean potential temperature. In theory, if R¢ is 
equal to unity then the buoyancy forces remove energy as fast as it 

is introduced by the shear, and the turbulence cannot maintain 

itself. However, this does not mean that R-¢ = 1 gives the criterion 


for the onset or disappearance of turbulence. 
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By substituting Eqs. (2.3) into (2.2) Rg, may be written 


as: 
kids 2.4 
Re = a Ri (2.4) 
where Ri is known as the gradient Richardson number and it is 
defined as: 
Rireses 90/ dz (225) 


8 (9dU/ dz)4 


Since the ratio Ku/Ky is not readily obtained from 
experimental data the gradient Richardson number (Ri), which can be 
identified as the ratio of buoyancy to inertia forces, is generally 
used in practice. It was shown by Batchelor (1953) that Ri is a 
relevant parameter for the specification of turbulent flow near the 
earth's surface under non-adiabatic conditions. The chief advantage 
of Ri is that it can be estimated from measurements of winds and 
temperatures at two levels. Observations have shown that Ri can be 
used as an index of stability but it has a disadvantage in that it 
varies with height and like R, it fails to predict the onset or 
disappearance of turbulence. Nevertheless, Ri has many practical 
applications and will be discussed further in connection with the 
log-linear law. 

Stability analyses indicated that there exists a critical 


Richardson number at which turbulence is damped out. However, the 
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correct value of this critical number is uncertain. Townsend (1958) 
and Kao (1959) found that the theoretical upper and lower limits for 
the critical flux Richardson number are 0.5 and -0.5 respectively. 
Laboratory measurements showed that in a stably stratified viscous 
boundary layer turbulence would not grow when Ri-> 0.0417 
(Schlichting, 1968). Furthermore, measurements in the atmosphere 
indicated that turbulence is not found when Ri> 0.2 (Lumley and 
Panofsky, 1964). Attempts were made by Rossby and Montgomery G935). 
Lettau (1949), Businger (1955) and others to replace Ri _ by a para- 
meter which would give a more complete description of turbulence. 


However, thus far, no entirely satisfactory parameter has been found. 
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CHAPTER III 
WIND PROFILE MODELS 


Sieh The Wind Profile Observations and the Power Law 

At the beginning of the twentieth century extensive experi- 
mental results for flow in smooth and rough pipes were obtained by 
Burgers, Zijnen, Hansen, and Nikuradse and summarized by Schlichting 
(1968). Based on these results it was found empirically that the 
velocity distribution in fully developed turbulent flow can be rep- 


resented by a power law in the form: 
Unicel. (330) 


where z is the distance from the boundary and p is an empirically 
determined exponent. Von Karman and Prandtl (Brunt 1934,>p. 236; 
Binder 1955, p. 122) found that p-=1/7 for Re< 5x 104 in 
smooth pipes and decreases to p = 1/10 for ~Re = 10° in rough 
pipes. 
Experiments to determine the velocity of the air close to 

a solid boundary began with Stanton (1911). Using a pitot tube he 
measured the air velocities to within 0.001 inches from the solid 
boundary and obtained a number of velocity profiles. In 1921 
von Karman analysed Stanton's observations and found that they fitted 
the one-seventh power law (p = 1/7) with great accuracy. 

The discovery that the velocity of the fluid flow in pipes 


can be represented by a power law prompted many meteorologists to 


lil ASTIAHD : 


21S00M SILTOAS GUIW - 


w81 tswod siz bas enoizevisedO sitter? balw oft oe 3 
~tisqxs sviensixs YruInso disltinews od to gninniged sd3 aA ” 
yd bemisido srew eeqig ree bas dioome at wold 10% esfueet [sjnem 
aniidrdoiidse yd bosiismmue bine. aha bis ,feecsH ‘ati 4S , e199 Tud 
sit asd3 ylisoixitgms bouot esw ti etivest eeed3 no beead . (8881) 
-qo1 sd ns0 wolt josludiui bsgo!lsvab vitut al qcottudizjalb v1isolsy 


‘ato? sf3 nt wel yswoq 8 yd beinses? 
([,€) Gs » 


yiisolzigqms ns ef q brs yusbruod si mort sonszaib afd at s sitsdw 
:0€£ .q , AECL aoutd) [sbaseyt bas namie soV | .dnsnogxs benimis3eb 
nt *O x 2 > 98 x02 ‘\l = q eds bavo® (SSI .q ,220l zsbale 
dguoz at Po = 9H .yot, Of\f = qv 02 esesaso9b bos sagig Azoome 

03 sealo zis sft to yitoolev sia snimis3sb of etmemiasqud - 
ed sdut 1031q & goieJ .(if@l) coans32 Atiw nagéd yrabruod bites & 
bifoe oi mot esdonit 100.0 midtiw ot setttoolev tis sd2 weit 
£S@l ni see! aboot to aden 5 boitedde a etibeavod 

bess? yortd Jers bauwd, baw enotseviesdo @ "none 32 bseyisas aS 7 


Dt asa 7 ay wr = 





tiie 

apply the power law to the measurements of the air flow near the 
earth's surface. As early as 1908 Akerblom studied the variation 

of wind speed with height using observations taken on the Eiffel 
Tower. Furthermore, he used the change in the wind direction between 
the top and the bottom of the Eiffel Tower to find the coefficient 

of viscosity of the atmosphere due to turbulence, 

The variation of wind speed with height became of increasing 
interest after the development of the theory of atmospheric turbulence 
by Taylor (1915a). Taylor (1915b) used pilot balloon observations, 
obtained over Salisbury Plain, England, in 1914, to calculate the 
skin friction of the wind blowing over the earth's surface. He ex- 


pressed the skin friction (F) on unit area of the earth's surface as: 


FWhe) ko 05 (3.2) 


where Q, is the wind velocity near the surface, p is the density 
of the air and k, is a constant which depends on the surface rough- 
ness and is determined from the wind profiles. For the ground at 
Salisbury Plain the value of kg, was found to be between 0.002 and 
0.003. 

The wind profile data were frequently required in the military 
operations, particularly in the chemical and biological warfare and 
in the problem of screening targets by smoke. Taylor (191 5c) Re witt 


the help of Cave, published numerous tables showing the wind velocities 
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at 1.25, 4 and 6 ft above ground at various stations. They wrote: 
"The observations described in the following 

pages were undertaken at the request of the Director 

of the Meteorological Office with a view to provid- 

ing information for the use of the officers responsi- 

ble for discharging of gas during a gas attack", 

According to Frost (1948), Hellmann in 1919 obtained wind 
profile data from recording anemometers set up at 5, 25, 50, 100 and 


200 cm above ground and fitted the observations to the power law in 


the form: 


U = Uy cok 3) 


where U is the mean wind velocity at height z and U, is the mean 
wind velocity at a constant reference height z, . Hellmann found 
that p = 0.27 when all of the observations were used and p = 0.33 
when only the night observations were used. 

Geiger (1927) suggested that the value of the exponent p 
is not constant but depends principally on height. He reasoned that 
because with increasing height the effect of ground friction diminishes 
the value of p becomes smaller. 

Heywood (1931) obtained wind velocities from two anemometers 
at heights of 12.7 and 94.5 m, above ground, over a period of two 
years. He found that the wind gradient between two heights close to 
the ground depends on the wind velocity of the upper of the two layers, 


the vertical temperature gradient, and the nature of the ground, 
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Fitting the observations to the power law of the form used by 
Hellmann, Heywood found that p = 0.26, which was in good agree- 
ment with Hellmann's value. 

Sutton (1932) emphasized that the value of p depends to 
a large extent on the vertical temperature gradient, a point which 
had previously been neglected. He used Heywood's observations and 
found that during the summer the change in p from midday to mid- 
night was from 0.07 to 0.17, while during the winter the correspond- 
ing change in p was from 0.08 to 0.13. Although, the change in 
p from summer to winter was slight it did show that p is least 
during the period of the strongest lapse and greatest during the 
period of the largest inversion. 

More conclusive results concerning the variation of p 
with stability were obtained by Barkat Ali (1932). He fitted the 
power law to wind measurements taken at Agra, India at heights of 
1.7 and 21.7 m above ground. The values of p were obtained 
separately for winter (Dec. - Jan.), hot season (Apr. - May) and 
monsoon season (Jul. - Aug.). It was found that during each season 
the value of p was largely determined by the existing lapse rate 
of temperature. In hot and monsoon seasons the value of p_ ranged 
from 0.0 with large positive lapse rates in the middle of the hot 
afternoon to 0.46 with extreme inversions at night. In the mid- 
winter season values of p as high as 0.87 occurred towards midnight. 

Further studies to determine the effect of lapse-rate of 


temperature on the value of p were carried out by Giblett (1935). 
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He grouped 50 ft and 150 ft wind speeds according to the temper- 
ature difference ( AT) between 143 ft and 4 ft. It was found that 
the ratio Uj,59/Usq increased with increasing stability from a 
value of 1.0 when AT =-5 F toa value of 1.17 when AT=OF 
and to a value of 1.98 when AT =/7 F. The power law was fitted to 
the U150/Uso ratios for each AT and the corresponding values 
of p were found to be 0.001, 0.143 and 0.630 respectively. 

Best (1935) studied the increase in wind velocity in relation 
to simultaneously occurring temperature gradients and found that 
the temperature gradient and the wind gradient are interrelated. 

In addition to dependence of p on temperature gradient, 
Paeschke (1937) found that p increases with increasing roughness 
of the surface. He carried out an experimental study and found that 
p varied from 0.20 over smooth snow surface with roughness length 
rae te Sycem, to)0,53 over se turniperield with oz = 45 cm. It was 
also found that p depends not only on the height, but also on 
the kind of plant cover. For example, over a wheatfield with 
z. = 130 cm, it was found that p = 0.29, while over a turnip field 
With 2. = 45 cm, Pea Oss. 

Carruthers (1943) summarized the relevant literature concern- 
ing the work that was done on variation of wind with height until 


that time. The summary of the result is shown in Table 3.1. 
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Table 3.1 


Values of p found by different investigators 


i 


Observation 
heights Value of p Source 
(m) 
Coe ee ee a a ae 
UME tov 3 OF LZ Fo.0733 G. Hellmann 1917, 1919 
Cae Om = 2 0.10 to 0.44 ALC best) 1955 
Various heights 0.14 £0.0.25 R. Geiger 1927 
Lad foe2 1. / O70 to 0.46 Barkat Ali 1932 
to 0.87 for large inversions 
Height not stated Biaalhe, Fou. ocarce 
Die CO ur oU OL J.S. Dines 
iSeeaor 9) C0ome co 0. L/ 0.G. Sutton (Heywood's data) 
1932 
ies > Wate Ate) 0.009- to 0.62 M.A. Giblett 1932 
Height not stated 0.009 to 0.62 0.G. Sutton 1934 


NS 


During the 1943-1945 period wind profile observations, up to 
a height of 400 ft, were taken at Cardington, England. The measure- 
ments were made at 5, 25, 50, 100, 200 and 400 ft above ground with 
instruments suspended from a balloon. Frost (1947) analysed these 
observations and showed that if the wind observations were grouped 
with respect to temperature difference between 400 ft and 5 ft then 
the mean winds fitted the power law and p was approximately a linear 


function of the temperature gradient. Further measurements of wind 
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speeds up to a height of 1000 ft were made at Cardington during 
1946. Analysing these data Frost (1948) considered only wind speed 
profiles for which the temperature difference between 4 ft and 1000 ft 
was between 5 F and 6 F and when the wind speed at any height between 
4 ft and 1000 ft was greater than 20 mph. He plotted logarithms of 
mean velocities of 64 observations at 8 heights against logarithms 
of heights and found that the points lie on a straight line having 
the slope equal to 0.149. Frost concluded that a value of p = 1/7 
is applicable to an atmosphere with a dry adiabatic lapse rate over 
the first 1000 ft, while p = 1/2 may be taken as appropriate for 
inversions. 

Johnson (1959) analysed wind profile observations obtained 
at Suffield Experimental Station, Alberta, and at O'Neill, Nebraska. 
The Suffield wind profiles, extending from 0.5 - 16 m, were measured 
over prairie grass and over snow surface while the O'Neill wind 
profiles were taken over short grass and extended from 2 to 16 m. 
He found that over the snow surface under near-adiabatic conditions 
pr= Oal65: Furthermore, it was found that the value of p increased 
with increased surface roughness and with increased stability. The 
maximum values of p for Suffield and for O'Neill were found to be 
0.75 and 0.4 respectively. 

A comprehensive study of wind speed profiles was completed 
by DeMarrais (1959). The data for this study were from the 125-m 
meteorological tower at Brookhaven National Laboratory. The study 


showed that the value of p decreased as the lapse rate changed 
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from marked stability to slight instability. Thereafter, the 

value of p leveled off under moderate superadiabatic conditions 

and then rose again in strongly superadiabatic conditions. It was 
also shown that the value of p was dependent on the wind speed 

and direction. The parameter p decreased with increased wind speed 
and it appeared that p was relatively smallest in the quadrant in 
which the wind had the longest unhindered fetch over the smoothest 
terrain. There was also an indication that the value of p decreas- 
ed with height. The uppermost layers had the smallest p-values during 
unstable conditions while the value of p was highest when the lapse 
rate was less than adiabatic. 

Panofsky, Blackadar and McVehil (1960) considered wind 
profiles in the height range of 11 - 46 m and produced a nomogram 
relating the power p to the roughness length Zo and to the inverse 
of the Monin-Obukhov length L. The nomogram, which was intended 
primarily for engineering applications, shows that the value of p 
increases with increasing values of Zo and decreases with increas- 
ing values of 1/L. 

Munn and Richards (1963) obtained the median values OL eED 
for each hour of the day from the 20-ft and 80-ft winds at Douglas 
Point, Canada, These p-values were plotted against time and the graph 
shows that the value of p was between 0.2 and 0.3 during the day 
and between 0.5 and 0.7 at night. 

Pettitt and Root (1965) examined wind measurements obtained 


from the 200-ft micrometeorological towers located at Montreal, Quebec, 
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and at Whiteshell, Manitoba. At Montreal, 35-ft and 200-ft wind 
speeds were considered in the analysis. The ranges of daytime p-values 
for winter, spring, and fall seasons were found to be 0.13 - 0.30, 
0.14 - 0.34 and 0.15 - 0.32 respectively, while the corresponding 
ranges for night-time p-values were 0.14 - 0.36, 0.21 - 0.35 and 
0.17 - 0.37. At Whiteshell 20-ft and 200-ft wind speeds were consider- 
ed and the values of p ranged from 0.18 - 0.20 in summer and from 
Opize= Un Oe lnewinter, 

Panofsky and Prasad (1965) showed that the power law exponent 


p given by: 


> = 2¢1n UD (3. 3a) 


oC 2) 


is auniversal function of z/z, and z/L’. where To etsethe 
gradient length and will be discussed in some detail in section 3.4, 
They found that the exponent p at a fixed height increases with 
roughness and with stability. Other studies were done to investigate 


the variations in p and will be discussed further in Chapter IV. 


Sack The Logarithmic Law 

While the power law representation of the variation of wind 
with height gave reasonably good results it had no theoretical foun- 
dation. The theoretical analyses of the phenomenon based on the 
equations of motion of fluids started with Prandtl in 1904. In his 
now classical paper presented to the Third International Congress 


of Mathematicians held in 1904 in Heidelberg Prandtl showed that for 
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a fluid of small viscosity, such as air or water, the viscosity will 
substantially affect the flow only in a thin layer adjacent to the 
surface. Prandtl called this layer the "Grenzschicht" (boundary 
layer). The boundary layer theory was extended after 1904 by Prandtl 
himself, von Karman, and others. By introducing a mixing length, which 
was assumed to be proportional to the height above ground, and using 
dimensional analysis, Prandtl was able to show that the velocity 
distribution in the turbulent boundary layer can be expressed by a 
logarithmic law. Von Karman (1954) wrote: ''The formulation of the 
logarithmic law was the end result of a long struggle to obtain 
correlation between theoretical ideas and experimental evidence". 
Although, the logarithmic law was primarily intended for 
studies of flow in pipes it was soon applied to the wind profile 
observations close to the ground. According to Frost (1948) , Chapman 
in 1919 found that the mean wind speed observations in the layer 
from 10 to 500 m were represented well by a logarithmic law in the 


form: 
U = A+B logz (3.4) 


where A and B are constants. In the same year Hellmann proposed 


another logarithmic law in the form: 
U = A+B log (2 +C) (3.5) 


where A, B and C are constants which vary with location and with 


time. 
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Rossby and Montgomery (1935) found that in an adiabatic 
atmosphere the wind velocity could be represented by the logarithmic 


law in the form: 


U = A [log (2 + 25) - log z, | (3.6) 


where A is aconstant and zg is a length which depends on the 
surface roughness. Over open grassland it was Found sthate 27) = op ence. 
Sutton (1953) gives the following forms of the logarithmic 


law based on Prandtl's theoretical work: 


uy, Z 


u 
For smooth flow: Utt= ~ In emit} + constant (3e7) 
Vv : 
“ Zz 3.8 
For fully rough flow: Ude=me—s In t= (3.8) 
k Zo 
Ux z-d 
For very rough surfaces: Ue rn in Gere Zane ee (3.9) 
fo) 
4 oe: ; 
where u, = ( T,/ ) is the friction velocity, 1, is the 


shearing stress at the surface, Z, is the roughness length, ¥V 
is the kinematic viscosity, k is von Karman's constant and d is 
the zero-plane displacement which is regarded as a datum level above 
which the normal turbulent exchange takes place. 

Thornthwaite and Halstead (1942) plotted the first and the 


second power of the wind velocity against the logarithm of height 
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and found that for night-time observations the first power of 
velocity was satisfactory to obtain a linear variation while for day- 
time observations the square of the velocity gave a better cit; 
Consequently, they formulated a wind law which was a combination of 


the logarithmic and the power law. This new law assumed the form: 
Ue =.) ((logeze-plog z.)) yy, loges yi/P (3.10) 


where log z, is the y-intercept and log a is the slope of a straight 
line obtained by plotting log z :U. The value of P was believed 
to vary between 2.0 with fully developed turbulence and some value 
less than 1.0 when turbulence reached its smallest actual value. 
Thornthwaite and Kaser (1943), using newly designed instruments, 
made precise wind profile measurements in’the 0.5 = 20 fC Layer. 
The data showed that the graph of Inz :U was a straight line 
during the morning and evening hours when the lapse rate was adiabatic. 
On the other hand, when the thermal structure was superadiabatic the 
graph was convex to the U-axis, and for inversions the curvature of 
the graph became concave. 
Halstead (1943) studied the Inz :U graphs and, by intro- 
ducing a stability term, extended the logarithmic law to fiteene 
wind profile observations made under all stability conditions. The 


modified logarithmic law, proposed by Halstead, has the form: 


Zz (elo h 
i, 28 + — (2) - 24) fecide 






















i - 
. ‘ 7 
Ss - 
- - - 
if * ; ~ - 
, a 7 
to tewog jeri? sia anorisvrsedo smit-sigin stot gsd7 bayed 
5 
-yab tot slinw noijsixsy 1seni! s stagdo oJ yrososteisse asw yotool 
a 
.3f2 sze3ted 6 svsg ytisolsv ofi Jo srsupe sid anobsevasedo amid 
70 soljsuidmos 8 esw doindw wel boiw s basalumi02d vss .VisnsupsemoD 
a 7 
‘mmo? sii bsmuees wel wen ate’ .wef tewoq a3 bas oimddiyegol fa * 
. 
\ 
aly q\ ly : 
(9 ; “> B gol | (,s gol - s gol)] = JU 
jdigisiya & Io sqola st : & 30f bas iqsstetni-y sii et ~s gol szadw 
7 S 
bs [od as G to suisv sAT J: s 30f gatttolq yd banissdo spit 
o emos bis sonsaludiut bagolaveb yilui dsiw 0. nsswised yisw oF 
.sulsv [sutos stesiisma eti bserossy sonsiudtu+s asdiw 0.1 gags ses! 
-8insmutjeni bengiesb yiwen goteu .(E#0]) sees. bob sitiswisnxodT 


| ° nN 
sI 33 8S - ¢.0 oft ni einsmstvaetsm slitorq baiw seissig sham 


snil ijidgisije © aBw U: s nf Yo"viqetg srit teat bewode sisb sit 
+ ” i *. 
- , 


2isedsibs esw stex-seqel oid oatlws2au0d ‘gnineys baa gninzom sid gataub 


= q » i 


e a> fe. > & oe . ae a ~ | : oe , } f 1 
383 2iIJBdSIDbsisqua 2sw stujoursse Jemaeds sis asdw ,basd xrsidoe ofa oP. 


a. 
— 


290 sIusgsvius Sd3 anolexsvnt sot bas aixs-U sit o3 xsvneo sew dqerg 


-9¥59n0> smsced tigatg 


-O1tni yd .bus adgets U: s mi arnt betbuse (EdOl) bsstaisH 






S 





+¢ 
; 


$43 312 o3 wel simdiisegol of4 bobnsixe ,mis2 ustitdsse s gatsub 
. - - a 


es > J 
Be a a ysilidste iis tebnu pp aR slite 
: 7 






imiod e. Bad , ‘beste isl yd beeogoxg 4s 
7 - wee ste 


' 7 : 
: A 


‘a 
Ss 








$ 


~256 


where Oy is the lapse rate of potential temperature, c is a 
constant and TS is the surface temperature. The stability term 

( Cap/To ) (2) - 21) is positive for stable air, negative for 
unstable air, and zero for neutral equilibrium. To test the valid- 
ity of the modified logarithmic law Halstead used part of the data 
obtained by Thornthwaite and Kaser and showed that the graph of 

In (25/21) : [ (Uy = lbp) = ca,/TQ) epee Z4)1 was, very 


nearly, a straight line for all stability conditions. 


3.3 The Deacon's Generalized Wind Profile Law 

A significant contribution to the formulation of the wind 
profile laws was made by Deacon (1949). Following an extensive 
investigation of the wind profiles in the lowest eight meters of the 
atmosphere Deacon proposed the generalized wind profile law inetne 


form: 


Ue ae (B12) 


where U is the horizontal wind speed at height z, and a and 86 
are constants in any given situation. The parameter 8B , whichere= 
flects the degree of curvature of the wind speed profile is a function 
of stability and assumes values: 8 > 1 for unstable conditions, 

8 = 1 at neutral stability and 8 < 1 for stable conditions. In 


the integrated form the Deacon's generalized wind profile law becomes: 


Uy z.1-8 
po | ee ab oie G.13) 
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For neutral conditions 8 = 1; and the equation (3.12) transforms 
into the logarithmic law. The parameter 8 , which is regarded as 
a measure of stability was evaluated by a number of investigators 
and several methods for its evaluation were developed. Deacon 
found that for velocity profiles over long grass the value of 8 


varied from 1.2 under superadiabatic conditions to 0.75 under marked 


inversions. The values of 8 were obtained from the equations: 


Z 
2a k= 6 
U Coe dueriel 
sal ee Pie ad ace ane (3.14) 
U Z, 1-8 
. (4) 0-1 
Oo 
and 
je ae Nes 
= wt =- -£ 
Uo ike wha Za 2 (3.15) 
Up = Uy 16 eas) 
29 3087 )4% 


where U,> U,; Uz are wind speeds at heights 2), Zo, 23 respectively 


and z is the roughness length. 


oO 
Rider (1954) used wind speeds at 37.5 cm and 150 cm, and a 
surface roughness of 0.32 cm, and found that the value of 8 ranged 
from 1.13 to 0.78. Furthermore, Rider discovered that the Deacon's 
generalized wind profile law was satisfactory under unstable conditions 
only and in order that the law would fit the observations under stable 


conditions, 2, would have to increase in magnitude with increasing 


stability. 
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Longley (1956) evaluated the parameter 8 using the 


relationship: 
U3 - U 
loges? ohates Z3 Z9 
B = Le Us Ce Uy 9 r = oe = ae (3,16) 
fozer ¢ 1 


where UL» Up» U3 are wind speeds at heights 2Z,, 29; Z3 respect- 
ively. The values of 8 were computed from several samples of data 
including data obtained at Suffield Experimental Station. A compar- 
ison of the results showed that it was impossible to obtain 8 
with sufficient accuracy to permit clear distinction between different 
types of stability and this led to the conclusion that its use in 
any practical problem was not warranted. 

Davidson and Barad (1956, 1957) analysed the data obtained 
at O'Neill, Nebraska and found that the parameter 8 varied with 
height. The variation was such that 8 decreased with height under 
stable conditions and increased with height under unstable conditions. 
They concluded that, " ... the Deacon formula, although adequate for 
gross wind interpolation purposes, does not accurately describe the 
effect of stability variation on the wind profile in the first six 
meters above the ground". 

Palmer (1956) calculated the values of 8 by using Barad, 
Craw, and Longley methods. The comparison of results showed that 
the values of 8, obtained by the Craw and Longley methods, were 


almost identical for the same heights. No direct comparison was 
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possible with Barad's method since it involved the use of wind data 
from different heights. It was concluded that the Longley method 
was superior to the Craw method since the former involved less 
computation and had less chance of an arithmetic error. 

Deacon (1956) suggested a reason for the failure of his 
generalized wind profile law to predict correctly the shearing stress 
under certain conditions. He found evidence that there was a marked 
rise in the drag coefficient as the wind speed decreased to a low 
value and this applied both to stable and unstable conditions. 
Consequently, there was an increase in the drag coefficient with de- 
creasing Reynolds number. Deacon concluded that the source of the 
reported failures of his generalized wind profile law was due to the 
absence of fully rough flow. 

In addition to evaluations of the parameter 8 several 
investigations were made to examine the parameters u,, 2, and d. 
Different methods for evaluating these parameters were given by 
Deacon (1949), Panofsky (1952, 1963), Kao (1959), Robinson (1962) 
and others. 

Deacon (1949) postulated that at sufficiently high fluid 
velocities the roughness parameter 2Z, is dependent only on the form 
and distribution of the roughness elements while at much lower fluid 
velocities the roughness elements are submerged in a shallow viscous 
sub-layer adjacent to the surface vets onl peat: v/ (9.05 uy) Deacon's 
observations over long grass showed a marked decrease of Zp, with 


increasing wind velocity and he concluded that the reason for this 
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decrease in z, was due to the long grass blades bending to the 

wind to produce a smooth surface. Based on Deacon's work Sutton (1953) 
prepared a table of z, and u, for natural surfaces. Sutton's 


table shows that Zo varies from 0.001 cm over mud flats to 9 cm 


over thick grass up to 50 cm high, while u, varies from 16 cm peone 
to 63 cm sec”! for the same two surfaces respectively. 

Hage (1961) yet iated zy from wind profile measurements 
obtained at Suffield Experimental Station and found that over sparse 
prairie grass Zz, increased with decreasing wind speed. The range 
of z, was from 2, = Gey) cneeateaUGuee= 6.45 m peck to z, = 7.4 cm 
at Up.5 = 1.00 m sec’! 

Dickson and Allbee (1967) evaluated the friction velocity 
u, for 1370 wind profiles using procedure outlined by Kao (1959). 

They found that the friction velocity is primarily affected by the 


velocity shear but not affected by thermal stratification. Thus they 


concluded that u, can be used in any atmospheric stability condition. 


3.4 The Log-Linear Law 

Most of the recent progress in the understanding of the 
variation of wind and temperature with height was brought about by 
Monin and Obukhov (1954). They introduced the, now well-known, 
similarity theory according to which there exist near the ground a 
velocity u,, a length L and a temperature T* that are essentially 
invariant with height. Under the conditions of similarity the wind 


profile in the surface boundary layer can be expressed as a function 


~@Qs 





etij o2 gnibnad esbaid asatg gnol sit o3 sub aw ot mt e2so798b 
(Eeer) mostue aaow e'moosed no boeed .soRiwwe dicome & ebybo1q 67 baiw 
a'nosju2 .e¢sontiwe I[swjen rot yu bus .s 206 sided & bexsqsxq 


m> ®@ ot at612 bum tovo mo 100.0 moz2 esizsv .s Jad3 swore sidad 


° 


[+ 


998 m3 Of movi estrav .uv ailtiw .fgid wo Of 03 qu sebzg A5ld? tsVvO 


.¥ievitosqesar nddadive ows emst of2 107 * 98 mo €8 63 

Sinema masem siiterq Sutw mors? >= Seiseleve (1ael) 93H 
sezsqe tevo jisdi bavol bne noitsss2 Iesnemiteqx® blaltive 38 bentaddo 
eget ott -b3eqe bniw gniasstosb diiw beeroezon) o* 88878 siilsexq 


zs io 


mo #.f = o* 09 358 m CA.d = 2.0" js moa \.0 = 3° mozt 2sw 6 


14508 m 00,i = 2.0" 38 
" etiey nofttols? oft batsulave (f8€l) eadJIA bas noedtoid = 
~(@C@L) o8% Yd barifiue siuvbs007%g caidd selhioww baiw OFEL fod 42 
eid yd badosiis yl isamizg ef yibeofav mottorst) sd3z sad? Baye? YoaT 


yeds aul .notteslitasiia Ismtsi? yo Se2z0F%s Jon gud eedde yitoolayv 


Moltibaos yititdete oFastigzomse \se 5k Beau sd nap yu taeda Bbebulonod 


7 
wet asgehi-gol siT ~.€ 

#43 20 galtbassexsbay a3 nd €#3igotg insset git to Jeo 
yd Juods ers enw tdap-ors Hj3iW swwosysqmes bos baiw io solisiqay 


atond-i lew won ,9f3 beoubortni yedT .(d20L) voraiudd bos nlm 


wt? eee jaixs sxad2 f2idw 02 gotbre228 yrosd3 adr 


r : ‘if 
lads | + on atnead aad f jatage 8 otzaley 
_ Pais ar4 ae ou ‘ bart om 

Ai ‘iii seals sie olan deal: ; 


> 


» 


’ 


OT; 
a 
4 

te 











230: 


of the non-dimensional height z/L in the form: 


S = ¢ (3917) 


poate ad 
S eee (3,18) 
and g = ¢(z/L) is a universal function. The Monin-Obukhov length 
L may be defined as: 
3 
pus Gy t 
aie cael (2219) 
k gi 


where T is air temperature, Cp is specific heat a constant pressure, 
H is vertical heat flux and the other symbols have their usual mean- 
ing. Once the function d is specified, Eq. (3.18) can be integrated. 

Considerable effort has been made over the last few years to 
find the universal function ¢ and several functional relationships 
were proposed. 


Monin and Obukhov (1954) suggested that for near-adiabatic 


conditions the function d can be expanded in power series as: 


d= ala seorpCeyr) +05 (2/L)* +£or* héal (3.20) 


ahs. 






:mtot efy nt JU\s #dated Ienotamemtb-non eda % 


~ 
-“ 
_ 
fm 

_/ 

a 

" 

co 


‘yd bentieb xsade boitw [anotenemtb-non 8 at @ szedw 


- 
(81.£) eS es Fe ) ae 
. Sob. LJ CJ - 
7 
i 
d3igasl vodaudO-ninoM sdT .nottonut Isexsvinu 6 at (I\s)& = a bos 
as beniteb ed ysm J — 
‘ee 
) 4U q ° 7. 7 
(@f.€) q 4 / 
- = — 
Hg a ' 


.Stvessig jnstenoo 8 json aitiosqe et g? ,31uIs79gmMS4 1is sit T srsdw - 
-1sem Isueu tisdt svsd elodmy2 isddo sd3 bos xuJi teed Isotiyev at H 7 ; 

-bojsigetnt od neo (81,€) .pd .bsitiosqe ei h mnolstonut eft son0 . gat 
OJ exs9y wet jeal sd3 revo sbem nesd esd srotts sidsrebtenod 


agidanolssiey {snoijonui Isisvse bas 4 nsoldonut Isexevinu sid batt 


. bs2ogozq eae i 


_ 


ees aan wot ted3 betesgegue (AZe!) voddudO bas inom om 


128 e9lise rewog nk bebusgxs od mso 8 notioqui oda enots “a 
¥ 


» 


— oe, | 
the Salas Seer i village 
— “4 a hee ; is rar 
, 7 7 ins ye 





=e 


then for small values of z/L the function can be approximated by: 
g = 1+4(2/L) (3.21) 


where @ is a constant to be determined from observations under near- 
neutral conditions. Substituting the function given by Eq. (3.21) 
into Eq. (3.17) and integrating yields the familiar Monin-Obukhov 


log-linear law: 


\ , fee, (3.22) 


Since the Monin-Obukhov length L can be measured only when 
the vertical heat flux is known, which is rarely the case, Panofsky, 
Blackadar and McVehil (1960) introduced a gradient length L' defined 
as: 


sa u, 9 (9U/98z) (3.23) 


kg (30/3z) 


where @ is the potential temperature and the other symbols have 


their usual meaning. The gradient length L' is related to, L. by: 
i a aL (3.24) 
where Ky and Ky are exchange coefficients for heat and momentum 


respectively. The question of the relative sizes of Ky and Ky 


has still not been answered satisfactorily. However, investigations 
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= 3D 
suggest that the ratio Ky /Ky is not too far from unity (Lumley 
and Panofsky 1964). 


Assuming that Ky/Ky ~ 1 the modified Monin-Obukhov function 


can be written as: 


g@ = I+ a(z/L') (3,25) 
and 
uy, 
ie ae I Sys ee are (3.26) 
k Ze Tee 


Another functional relationship which was found particularly 
useful for specifying the function d was derived in different ways 
by Kazansky and Monin (1956), Ellison (1957), Yamamoto (1959), Panofsky 


(1961), and Sellers (1962), and has the form: 


we = 1B (2) Ho) sil (Sp27) 


This relationship became known as the KEYPS function after the initials 
of its inventors. 

Priestley (1960) advanced the view that the values of the 
Monin-Obukhov coefficients (a's) are products of smoothing rather than 
true physical constants of a steady regime intermediate between free 
and forced convection. He defined a junction height between free and 


forced convection as: 


getey) 3 
$ = = ke bh (Ky /Kyp) ! (3.28) 
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BIE, 
where for k=0.4, h=0.9 and Ky/Ky = 1, the junction height 
z/L = -0.0316 at the unsmoothed junction, Furthermore, Priestley 


specified the function ¢ for free convection as: 


Knee : 
g- = dE) 7 Be (3.29) 
M 


where h is the Priestley constant. 
Webb (1960) suggested a two-part function ¢ with continuous 
first, second, and third derivatives at the junction. This function 


tsporethe form: 
Gaee lets 5(2/ UL") forme- 2/1) eoeOsO sym ( o.0) 
d = 0.316(-2/t')"2/3- 0.00143(-2/L") "4/9 for -2/L' > 0.0317 . (3.31) 


Businger (1961) studied the spectrum of atmospheric turbulence 


and proposed a universal function in the form: 


K 
48/3, ee 18/3. Mbeos tetas hy Ne (3.32) 
Cc 


where K,/K, = 1.7 and Rg is the flux Richardson number. 
Panofsky (1963) showed that, in general, the log-linear law 


can be written in the form: 
u 
k 


U = [In > - W(z/L') ] (3433) 
te) 
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where v (z/L') is another universal function which is related to 


d(z/L') by: 


v (2/L') or es ee dé (3. 34) 


Brooks et al (1963), following Sheppard's (1958) arguments, 
found that a better fit to observed wind profiles was possible by 
applying a square root modification to the Monin-Obukhov universal 


function and defined a new universal function as: 
z'\% 
a ae | (3.35) 


where’ 2)’ =2-d and-z' >z.. In the integrated form the Brooks 
log-power law becomes: 
Us, eA 


U2") = = [In 2 + 27, @) ] (3.36) 
@) 


where Yj was found to be 1.24. 

Swinbank (1964) introduced a new non-linear height variable 
X which had the property that the graph of U: log X was a straight 
line no matter what the stability. He then derived a universal function 


in the form: 


pux = [1 - exp(-Z/L) ee ‘oe 
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which upon integration yields the Swinbank exponential wind profile: 


exp (Zo /L) - l 


u 
* 
Uy - U; rd Si ae ln Spa I (3.38) 


where the symbols have their usual meaning. 
McVehil (Lumley and Panofsky 1964) found that in stable air 
for Z/L'< 0.3, (Ri < 0.1), the log-linear law assumes the form: 


tives GZ ea 9 
cael ey (3.39) 


where B= 1/Ri - L'/Z. However, when Z/L'> 0.3 no simple 
relationship exists for the universal function p . Furthermore, it 
was suggested that at large values of Ri turbulence changes character 
and is gradually replaced by gravity waves. 

Brooks, Mulholland and Smith (1965) experimented with a 
hyperbolic and an inverted exponential form of the universal function. 
This resulted in a log-law plus power law binomial given by: 


u(zZ') = Ae TEN (3.40) 


where Z, is a characteristic height, b is the exponent and c 
is a constant. 
Pandolfo (1966) established the form of the universal function 


in the constant flux boundary layer in lapse stratification. He found 
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that Monin-Obukhov function applies in lapse forced convection 
conditions and Priestley's function applies in the free convection 
conditions. The integrated wind profile laws for forced and free 


convection then become: 


u Z, 
U(Z9) - U(Z{) = - fn - + 2 (2) - 21) ] for 0 > > @), 
| (3.41) 
uy I, 7, : : 
U(Z>) = U(Z,) = = 6 &)? i I—|1/6_ | — 1/6) for (Z/L) , SSA nt 
(3.42) 


4 -2 , ; 
where c = 3k /3 h /3 and h is the Priestley constant. 

With the introduction of the log-linear law it became necessary 
to determine the constants which enter into the equation. In addition 


to friction velocity u, and roughness length z which entered 


°°? 
into the logarithmic law, the log-linear law requires the Monin-Obukhov 
constant oa andthe length L. In the Monin-Obukhov original analysis 
it was found that a= 0.6. However, it was shown by a later study 
(Taylor 1960) that the value of © is about six. The reason for this 
discrepancy was ae to the fact that Monin and Obukhov applied the log- 
linear law to a larger range of stability than was justified. 

Taylor (1960) calculated the values of ao from the data 
given by Monin and Obukhov, Swinbank, and Rider and obtained different 
values of © for forced convection, free convection and inversion 


conditions. The numerical results obtained by Taylor are given in 


Table 3.2. He concluded that the value of a is about six for 
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stability corresponding, approximately,, with forced convection 


conditions. 


Table 3.2 


Values of a 





Monin- 
Rider Swinbank Obukhov 
Forced convection (Lay eo Ot 6,1) fal 6.072 039 
Free convection OFS52)0,035 1,432 0.14 e314 5 0505 


Inversion De wt) e5 get mals) -- 


ee ______ EE e 


Priestley (1960) found that the values of the Monin-Obukhov 


constants in superadiabatic conditions can be predicted by the equation: 


: : i 
Oj; = (eh ya (0.7121) (3.43) 


(Si 


where y= Ky /Kyp> k = von Karman's constant, h = Priestley's 
constant and i= 1, 2, 3,... . The numerical values of 4, then 
become a = 5.63, % = 36.3, Og = 191, and so on. When the 
coefficients given by Eq. (3.43) are substituted into Eq. (3.20) the 
function d consists of terms of alternating sign since L is 
negative for superadiabatic conditions. Consequently, the function 


¢@ converges for all finite 2. 
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McVehil (1964) analysed wind and temperature profiles from 
O'Neill, Nebraska, and Antarctica and found that in stable air the 


value of a is approximately seven. 
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CHAPTER IV 


ANALYSES AND RESULTS 


4.1 The Power Law 

The current analyses began with the determination of the 
power law exponent p. By differentiating the power law given by 
Eq. (3.3) logarithmically it can be shown that, in the finite differ- 
ence form, the power law exponent is given by: 


~ Kins 
oi iea Cy 


By plotting in z: ln U and using the least squares method the value 
of p was calculated for each profile. In addition, to determine the 
accuracy with which the observed data are represented by the power 

law the correlation coefficient between 1nz and I1nU was deter- 
mined for each profile. The logarithm of height was considered as the 
independent variable. The application of eeeee reece ton coefficient 


will be discussed further in the next chapter. 


In order to compare the values of p obtained in this study 
with the values obtained by different investigators the method of 
comparison used by DeMarrais (1959) was adopted. The wind profiles 
were grouped into the following categories: superadiabatic, adiabatic, 
stable and inversion. The average value of p was then obtained for 
each category. The values of p obtained in this investigation are 


compared with those obtained by other investigators in Table 4.1. 
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Table 4.1 


A comparison of the values of p obtained in this study with those 
obtained at other locations 





Height 
Site Range Terrain Super- Neutral Stable Inversion 
(m) adiab. 
Suffield 0.5-16 OFZ 2 0.23 -- 0.31 
(41 profiles) 
Suffield 2 - 92] sparse prairie 0.18 Oke ‘Ohyalv Wl EP 
(129 profiles) grass 
Suffield 16 - 92 0209 0.14 O79 0.36 
(129 profiles) 
Kerang, 0.5-16 grazing land Oy HE -- -- -- 
Australia 
Suffield(1959) 0.5-16 prairie grass -- -- Oez2 0.49 
Suffield(1959) 0.5-16 snow 0.16 Oel7 -- 0.28 
O'Neill, Neb. 2 - 16 short grass Om -- QO. 17 0.38 
Quickborn, 10 - 70 meadows 0325 One -- 0.61 
Germany 
Tallmadege, Ti- 49 eetlatertield 0.16 0.20 05.25 0.36 
Ohio 
Hanford, 15 -122 mountainous 0.09 OF 12 0.14 O25 
Washing. 
Cardington, 8 -120 grass field 0.15 ORL, O74) OR oZ 
England towO. a, 
Harwell, 9 - 27 airfield 0.09 0.08 0.18 -- 
England 
Idaho Falls, 6 - 61 #£4desert 0.15 0.18 Op22 -- 
Idaho 
Brookhaven 11 -124 | nearby OnLo 0,29 Os39 0.46 
COOso2 
Brookhaven 11 -125 |wooded OneLa O728 0738 O72 
Brookhaven 46 -125 area 0.18 Oss) 0.45 Oper 
Brookhaven 11 - 46 ae |» takes 0.34 0.37 


i ereeee en 
* Data obtained in the present study 
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A study of the results in Table 4.1 reveals that the value of 
p depends on the location, stability, height range, and surface rough- 
ness. It is observed that the values of p obtained for Suffield 
0.5-16 m layer by Johnson (1959) are in fairly good agreement with 
those obtained in the present study. The small differences which do 
exist may arise from the fact that Johnson calculated the values of 
p for a specific temperature gradient while in the present study 
average values of p are given. In the Suffield 2-92 m layer the 
values of p appeared to increase in extreme superadiabatic conditions. 
Consequently, a relatively high average value of p was obtained for 
lapse conditions. Such rise in the values of p under extreme super- 
adiabatic conditions was also reported by DeMarrais (1959). 

Investigations of wind profiles carried out at the Brookhaven 
National Laboratory indicated that, in addition to the temperature 
gradient, height range and surface roughness, the value of p 
varies with the vertical wind shear. In order to account for the 
variation of wind with height and to obtain a measure of convective 
stability the gradient Richardson number was calculated for each 
profile. The calculations were performed using the gradient Richardson 


number in the finite difference form; defined as; 
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The potential temperature 9 was calculated from Poisson's equation 


in the form: 


/ 
1000 8 
@ = T Goa. a) (4.3) 
where T is the air temperature, R/C, = 0.286 for dry air and 


P 


P is the barometric pressure. 


Next, wind profiles in each of the four height ranges shown 
in the first four lines, in Table 4.1 were grouped in decreasing 
order with respect to Ri and divided into nine subgroups in such 
a way that there were approximately equal numbers of wind profiles 
in each subgroup. The average value of p and the average value of 
Ri were then calculated for each subgroup. Graphs of p : Ri 
were constructed for each height range and are given in Figs. 1-4. 


The discussion of the graphs will be given in Chapter VI. 
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BIG. AL. Variation of the power-law index p with 


Richardson number in the 0.5 - 16 m layer (Kerang 
data). 
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Big. 2, Variation of the power-law index p with 
Richardson number in the 0.5 - 16m layer (Suffield 
data). 
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FIG 3’. Variation of the power-law index p with 
Richardson number in the 2 - 92 m layer (Suffield data). 
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FIG. 4. Variation of the power-law index p with 


Richardson number in the 16 - 92 m layer (Suffield 
data). 
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4.2 The logarithmic Law 
The logarithmic law in the form given by Eq. (3.8) was 
considered as appropriate for the sites at Suffield and at Kerang. 
For the Kerang data the parameters uy and 2Z, appearing 
in the logarithmic law were evaluated by Swinbank (1964). Hence 
it was necessary to evaluate uy and Zo only for the Suffield 
terrain. The method of evaluation was similar to that outlined by 
Panofsky (1963). By taking the derivative of the logarithmic law 
in Eq. (3.8) it can be shown that in the finite difference form 


u, is given as; 





wy = ky (4.4) 


Furthermore, by applying the boundary condition U=0 at z= zy 


it can be shown that z, is given as the z-intercept of the re- 
gression line obtained by plotting 1nz: U. Graphs of Suffield 
wind profiles from which u, and Zz were calculated are shown in 
Figs. 25-30, Appendix A. 

Once the values of z, were obtained from graphs the 


values of u, were recalculated using the drag coefficient equation: 


* 


Ree wo k (4.5) 


-2p- 
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The value of von Karman's constant k was assumed to be 0.4 which is 
the presently accepted value. The values of u, and Zn calculated 


from the Suffield wind profiles under neutral or near neutral conditions 


are given in Table 4.2. 


Table 4.2 


Calculated values of u, and z, for Suffield 





a Th"T9 5 Uo.5 Zi u, (m sec’) uy, (m sec”) 
profile (F) (m sec” 1) (cm) (graphical) (drag coeff.) 
12/2/60 -0.4 Bal 9 0.48 0.56 
13/1/64 -0.3 hate 323 Om25 0.25 
23/7/03 -0.3 2710 3.4 0.38 0.36 

1/2/60 -0.3 2.8 2.8 Dees 0.38 
14/10/65 0.0 2.8 6 0.48 0.38 

8/2/60 0.0 Bhal 2.9 0.45 0.46 


a 


The exact nature of the surface at the times when the profiles 
were obtained was not available. However, the date of each profile 
gives some indication of the surface roughness, discussed in Chapter I. 
It can be observed from Table 4.2 that there is considerable scatter 
in the calculated values of uy, and Z,- By assigning more weight 
to the profiles for which the correlation coefficient between ln z 


and U was higher the average values of u, and Zo for Suffield 


terrain were found to be 0.34 m Spr and 3 cm respectively. 
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It was pointed out by Lettau and Davidson (1957) that a 
vegetation cover of varying density can raise the reference level at 
which the boundary condition U=0 is satisfied significantly above 


the average physical soil surface. Under such surface conditions z 


fe) 
can be determined more accurately by using Eq. (3.9) in the form: 
Uy 
Zz + D 
U = i in Hy (4.6) 


oO 


where D = d+z, is the zero point displacement parameter. The 
boundary condition now becomes U=0 for z = -d. It was shown 
by Robinson (1962) that, given at least three observations of U 


anumezo mld... (4.0) can be solved? for “uy, .7z and D using a 


oO? 
secant iteration method. A computer program prepared by Robinson 
was obtained from the Department of Meteorology University of 


Wisconsin, and values of u,, z and D were calculated using 


o? 
the Suffield data. The results of these calculations are presented 
in Table 4.3. 

Comparison between Tables 4.2 and 4.3 shows that, in general, 
the values of u, are in agreement with one another. On the other 
hand, the values of 5 in Table 4.3 are, approximately, only half 


as large as those in Table 4.2. The values oft Sul ys and D 


Oo? 
obtained from the wind profile dated 14/10/65 are somewhat puzzling, 


particularly the positive D value. Such anomalous D values were 


also reported by Lettau and Davidson (1957) when they analysed the 
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Great Plains Turbulence Field Program Data. By omitting the anomalous 


wind profile from Table 4.3 the average values of u,, 2 and D 


oO? 


were found to be 0.34 m sec", 1.5 cm and -3.2 cm respectively. 


Table 4,3 


Calculated values of uy, Z,, and D for Suffield 


i 


vee Ty-To 5 a Ze D 
profile (F) (m sec” |) (cm) (cm) 
12/2/60 -0.4 0.44 1.0 -3.7 
13/1/64 -0.3 0.22 VG -1.2 
23/7/63 -0.3 0.35 hae -5.4 

1/2/60 -0.3 081 1.7 263 
14/10/65 0.0 0.65 19.4 64.5 

8/2/60 0.0 0.40 135 =o75 


nc nn eT 


4.3 The log-Linear Law 
as The Monin-Obukhov log-linear law 

The Monin-Obukhov log-linear law given by Eq. (3.26) was 
considered in the analyses. Because heat flux measurements were not 
available the gradient length L', given by Eq. (3.23), was calculated 
for each profile. The Monin-Obukhov constant a was calculated from 
the wind profiles in each height range by plotting the quantity 
In z + d(z/L') as a function of wind speed and using the least squares 


method. It was found that the value of o was different for free 


~9a~ a 
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Oe 
convection, forced convection, and inversion conditions. Moreover, 
the value of 4% was different in different layers. The calculated 


values of a for each height range are given in Table 4.4. 


Table 4.4 


Calculated values of a 





Height range a Range of Ri 


Suffield 2 - 92m 0.058 Ri < -0.500 
0.58 =O, 00a Ria 0 
4.5 ORseRie<n0 eLa0 
ae) Ri > 0.140 
Suffield 16 - 92 m 0.058 Ri < 0 
1.,.0 0«<-Ri—<-0.2390 
0.058 Riv 30.390 
Suffield 0.5-16 m 0.058 Rice 
6.0 Ries 
Kerang 0.5-16 m 0.58 Rice) 


EEE 


The values of @ in Table 4.4 bear some resemblance to those 
obtained by Taylor (1960) for forced convection, free convection, and 
inversion regimes. Moreover, the value of % obtained from Kerang 
data is in agreement with the value a = 0.6 used by Bernstein (1966) 
when he used the Kerang data to compare different forms of the log- 


linear law. 
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b. The Brooks log-power law 

Following the Brooks et al (1963) suggestion that a better 
fit to observed wind profiles was possible by using a square root 
modification of the log-linear law, the Suffield and Kerang data 
were fitted to Eq. (3.36). As in the case of the Monin-Obukhov log- 
linear law the gradient length L' was used in the analyses. The 
parameter Y' which appears in the Brooks log-power law was cal- 
culated from the wind profiles in each height range using the least 


squares method, The calculated values of y' are given in Table 4.5. 


Table 4.5 


Calculated values of Y' 


EEE 


Height range ye Range of Ri 
ee a a cl ee een en eens at oe 
Suffield 2a 925m 0 R70 
0.1 Ri < 0 
Suffield 16 - 92m '.5 Ri > 0 
Oo Rivs20 
Suffield 0.5-16 m li. 5 Riu 20 
Q,.1 Ri < 0 
Kerang 0.5-15 m la, O Ries 0 


en EEE Sean 


A study of Table 4.5 reveals that the parameter y' plays 
a similar role in the Brooks log-power law as the parameter 4 
does in the Monin-Obukhov log-linear law. However, it appears that 


Y' does not vary with the thickness of the layer considered. 
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Cc. The modified log-linear law 

After some experimentation with different wind profile 
models discussed in Chapter III it was found that the data were 
represented well by a modified log-linear law in the form: 

deecietele (15h 2a eucbe Shite) (4.7) 
k Zo te JL 

Here L, is a constant length determined from the wind profile data 
and L is the gradient length given by Eq. (3.23). Calculations 
showed that Lo is relatively small for Ri >O and large for Ri <0. 
Furthermore, Ly is different in each height range. The calculated 


values of ys for each height range are given in Table 4.6. 


Table 4.6 


Calculated values of L 


nena EEE SARNEREnE 


Height range Le (m) Range of Ri 
eR 
Suffield 2 - 92m Sez > RiceeO 
Sy 2ax LO Ri < 0 
Sutrield: \L6.= 92°m 25 9 Rive: 0 
25° xKeL0 Ricea0 
Suffield 0.5-l6m 0.07 3 Ricoao 
TS .bex aL Rigas0 
Kerang 0.5-16 m 30 Rivcs0 
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It may be observed from Table 4.6 that there is a relatively 
large difference in US between stable and unstable conditions. 
However, it should be noted that in Eq. (4.7) the numerical value of 
(1/,)7 is of importance and not L, itself. 

The question now remains which of the wind profile laws fits 


the data most accurately. This question will be discussed in the 


next chapter. 
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CHAPTER V 


COMPARISON BETWEEN THE DIFFERENT WIND PROFILE MODELS 


One of the fundamental questions in micrometeorology which 
still awaits the answer is to determine which of the wind profile 
models describes the actual conditions best. Relatively few compar- 
ative studies between different wind models were done in the past. 
The method of comparison, generally, varied with each investigator 
and until the present time there is no standard way of making the 
comparison. 

Frost (1948) compared the power law with the logarithmic law. 
The results showed that over open grassland, for heights below 10 m, 
the power law representation was better than that given by the log- 
arithmic law. The wind-speed values given by the logarithmic law in 
the layer below 10 m were too low. In the 10-100 m layer the wind 
speeds were represented equally well by the two laws. On the other 
hand, in the 100-800 m layer the logarithmic law was better while the 
power law gave too high speeds, When the roughness of the surface 
was taken into account, it was found that for relatively smooth 
surfaces, with 2, « l cm, the wind profiles can be best represent- 
ed, up to several hundred meters, by a logarithmic law. However, for 
very rough surfaces the power law gave a better representation. 

Johnson (1959) analysed wind profiles extending from 0.5-16 m. 
The measurements were taken over prairie grass and over snow surface. 


It was found that under all stability conditions and for each surface 
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roughness the power law gave equal, and in most cases better, results 
than the logarithmic law. Furthermore, the results indicated that 
under all stability conditions the wind speed increased with height 
more rapidly than predicted by the logarithmic Law. 


DeMarrais (1959) considered the power law in the form: 
U = U, (—)P 
= Uy, ) @.1) 


He then used the observed wind speed at 11 m and the value of p 
for each ten-minute period to compute the wind speed at 46 m. Next, 


he computed the 46 m wind speed using the logarithmic law in the form: 


log (z3/z) 


Uy ee (U, - U)) Tog (23/21) (5 32) 


where the subscripts 1, 2, and 3 refer to 11, 46 and 125 m respectively. 
The computed 46-m wind speeds, for each law, were then compared with 
the observed wind speeds and the mean error for each temperature lapse- 
rate class was computed. The results showed that the logarithmic law 
had greater error than the power law for all except the superadiabatic 
lapse rates. 

Volkovitskaya and Mashkova (1963) analysed the wind data 
obtained from the Russian 300-m tower. They found that in the bottom 
25-30 m layer the wind speed varies according to a logarithmic law 


while in the overlying 25-300 m layer the wind profile deviates from 
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the logarithmic law and can be approximated by a power law. 

Thuillier and Lappe (1964) analysed the data obtained from 
the 1420-ft Cedar Hill tower located near Dallas, Texas. They found 
that for near-adiabatic and slightly unstable conditions the logarithmic 
law represents the data well up to a height of 300-400 ft. Above this 
height the wind speed was found to be nearly constant with height. 
For stable non-inversion conditions the power law was found to fit 
the data better than the logarithmic law while under inversion 
conditions the wind profiles were highly variable and dependent on 


the nature of the inversion present. 


The main area of research during the recent years involved 
the comparisons of different forms of the log-linear law. Panofsky, 
Blackadar and McVehil (1960) plotted the non-dimensional shear S 
versus ln WZ | Us where y was assumed to be a constant, and compared 
the theoretical relationship between S and yz/L, for KEYPS and 
Monin-Obukhov functions, with observed values obtained at five 
different locations. The comparison showed that the KEYPS function 
fits the observations in unstable air better than the Monin-Obukhov 
function. However, on the stable side neither function fits the 


observed data well. 


Panofsky (1963) made a theoretical comparison between the 
Monin-Obukhov, Webb, and KEYPS profiles by plotting the universal 
function ln y: ln z/L. The comparison showed that for small values 


of Inz/L_ the three profiles are identical. However, for 
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In z/L> 0.1 the KEYPS and Webb profiles agree well with each other 
while the Monin-Obukhov profile diverges considerably. 

Bernstein (1966) used the data obtained at O'Neill, Nebraska, 
and at Kerang, Australia, and compared the Monin-Obukhov, KEYPS, and 
Swinbank profiles with actual observations. The comparison was made 
by plotting the logarithm of the non-dimensional wind shear S&S 
against the logarithm of -2/L under unstable conditions. He found 
that the data were not of sufficient accuracy to distinguish among 
the different hypotheses and concluded that the presently available 
data are not sufficient to verify or refute any of the three models. 

Hansen (1966) compared nine wind-profile models for the 
diabatic boundary layer (including Monin-Obukhov, Brooks, KEYPS, Webb, 
Businger, and Swinbank models). The data used in the analyses were 
from O'Neill, Nebraska, Kerang, Australia, and White Sands Missile 
Range, New Mexico. The method of comparison was similar to that of 
Bernstein (1966). Hansen found that in the forced convection regime, 
0 > Ri > -0.02, all models allowed accurate determination of the 
roughness length and would give good estimates of the stress and 
vertical heat flux. However, in the transition zone between forced 
and free convection none of the models gave good results beyond the 
junction height Ri = -0.0317. Furthermore, in the free convection 
regime, Ri < -0.05, mone of the models fitted the experimental data. 
For the White Sands data, which was obtained under non-stationary 
conditions owing to the non-uniform terrain and the formation of the 


internal momentum boundaries caused by roughness discontinuities, 
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none of the models was found to be valid. Consequently, Hansen 
concluded that the present surface boundary layer models have limit- 
ed use in less than ideal situations. 

The comparisons between the power law and the log-linear law 
are rare. The main difficulty is that the former law is empirical 
while the latter is theoretical. This results in the absence of common 
factors, making a direct comparison difficult. 

The method of comparison adopted in this study attempts to 
determine the accuracy with which the least-squares profiles represent 
the observed data. This method is somewhat similar to that used by 
McVehil (1964). For each of the four height ranges shown in Table 4.4 
the wind profiles were grouped in decreasing order according to the 
gradient Richardson number into eight or nine subgroups, depending on 
the number of samples. The grouping was done in such a way that there 
were approximately equal numbers of samples in each subgroup. The 
average value of Ri was then obtained for each subgroup. Next, 
for each subgroup the correlation coefficients, discussed in the 
previous chapter, were transformed into Fisher's z' - statistic 
and weighted before being combined. The equations used in transform- 
ing the correlation coefficients are given in Appendix B. A procedure 
similar to that outlined by Brooks and Carruthers (1953) was used to 
apply the weights. This procedure was repeated for each wind-profile 
law. The results of comparison between the power law, the logarithmic 
law, and the Monin-Obukhov log-linear law are summarized in Figs. 5-8. 


Similar method of comparison was used to compare the three log-linear 
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wind-profile models discussed in Chapter IV, section 3, The results 
from this comparison are displayed graphically in Figs. 9-12. These 
graphs will be discussed further in the next chapter. 

The next question of importance in the analyses was to 
determine if the differences between different wind profile laws 
were significant. The five per cent level confidence intervals 
were considered appropriate to determine the significance, The 
significance tests were performed in two parts. 

First, for each subgroup, in addition to the weighted corre- 
lation coefficient, the five per cent level confidence intervals 
were calculated. The equations used in calculations are given in 
Appendix B. Figs. 13-16 show the comparison between the five per 
cent level confidence intervals for power law, logarithmic law, and 
Monin-Obukhov log-linear law. In addition, the comparison between 
the confidence intervals for Brooks log-power law, modified log- 
linear law, and Monin-Obukhov log-linear law is given in Figs. 17-20. 

Second, instead of computing the weighted correlation coef- 
ficients and the corresponding confidence intervals for each subgroup, 
an average weighted correlation coefficient together with the five 
per cent level confidence intervals was calculated for each of the 
four sets of data. The comparison of significant differences between 
power law, logarithmic law, and Monin-Obukhov log-linear law, for 
each set of data, is displayed graphically in Figs. 21-24 while the 
comparison between Brooks log-power law, modified log-linear law, 


and Monin-Obukhov log-linear law is given in Table 5.1. 
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Table 5.1 


Comparison of significant differences between Brooks log-power law, 
modified log-linear law, and Monin-Obukhov log-linear law 





Wind 
profile 
law 


Brooks 
log- power 
law 


Modified 
log-linear 
law 


Monin- 

Obukhov 

log-linear 
law 


ne LEE EES SSS 


Height 
range 


(m) 


Suffield 2-92 
Suffield 16-92 
Suffield 0.5-16 


Kerang 0.5-16 


Suffield 2-92 
Suffield 16-92 
Suffield 0.5-16 


Kerang 0.5-16 


Suffield 2-92 
Suffield 16-92 
Suffield 0.5-16 


Kerang 0.5-16 


No. 
of 
obs. 


Average 
weighted 

correlation 
coefficient 


0.985 


0.985 


0.990 


05999 


0.988 


0.985 


0.996 


05999 


0.987 


0.987 


05993 


e999 


5% level 
confidence 
intervals 


0.982-0.988 
0'.981-0.959 
0.985-0.994 


0: 998-0..999 


0.985-0.990 
0.981-0.989 
0.994-0.998 


0,.998-0.999 


0.984-0.989 
0.982-0.990 
0.989-0.995 


0.998-0.999 


This concludes the analyses of the wind profile observations, 


and in the following chapter the assessment of the results, here 


presented, will be made. 
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FIG. 8; Comparison between Monin-Obukhov log-linear 
law, logarithmic law,and power law in the 16 - 92m 
layer (Suffield data). °Ri = -3.084, [Ri = -11.479 
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FIG, 10; Comparison between Monin-Obukhov log-linear 
law, Brooks log-power law, and modified log-linear law 
in the 0.5 - 16 m layer (Suffield data). °Ri =-0.129 
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law, Brooks log-power law,and modified log-linear law in 


the 16 - 92 m layer (Suffield data). °Ri =-3.084, 
ORi = -11.479 
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BUG. (13. Comparison between the five per cent level 
confidence intervals for different wind-profile laws in 
the 0.5 - 16 m layer (Kerang data). 
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FIG. 14. Comparison between the five per cent level 
confidence intervals for different wind-profile laws in 
the 0.5 - 16 m layer (Suffield data). °Ri = -0.129 
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the 2 - 92 m layer (Suffield data). °Ri = -0.960 
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FIG. 18. Comparison between the-five per cent level 
confidence intervals for different wind-profile laws in 
the 0.5 - 16 m layer (Suffield data). °Ri = -0.129 
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Fig. 21. Diagram to illustrate the significant diff- 
erences between different wind-profile laws. The shaded 


areas represent the five per cent confidence intervals 
for the 0.5 - 16 m layer (Kerang data). 
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Bie. 22. Diagram to illustrate the significant diff- 
erences between different wind-profile laws. The shaded 
areas represent the five per cent confidence intervals 
for the 0.5 - 16 m layer (Suffield data). 
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Figs «23. Diagram to illustrate the significant diff- 
erences between different wind-profile laws. The shaded 


areas represent the five per cent confidence intervals 
for the 2 - 92 m layer (Suffield data). 


1.000 


0.990 
seesssderessnargaaccdsssreeeeeessessssoeesesseaseeesscees 
pepe 


0.980 


ORE WAS 54: Power law 


Logarithmic law 


0.960 Monin-Obukhov log- 
linear law 





0.950 
-1.00 -0.50 0.00 0.50 1.00 iPesy 6 


Ri 


Fig. 24. Diagram to illustrate the significant diff- 
erences between different wind-profile laws. The shaded 
areas represent the five per cent confidence intervals 
for the 16 - 92 m layer (Suffield data). 
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CHAPTER VI 
DISCUSSION OF THE RESULTS 


Before discussing the results it is desirable to consider 
some of the obvious sources of error. The Suffield wind-speed 
measurements, as well as temperature measurements, were given to 
one decimal place. It was estimated that at Suffield winds and 
temperatures were measured with an accuracy of about 5 cm Bec and 
0.05 F respectively. According to Swinbank (1964), winds and temper- 
atures at Kerang, which were given to two decimal places, were measured 
with an accuracy of about 2 cm sec~! and 0.02 C respectively. Further- 
more, as pointed out by Hage (1961), the Suffield wind- speed measure- 
ments were made with cup anemometers and are likely to contain 
systematic errors because of turbulence. In addition, in the early 
16-m tower experiments, the top anemometer (16 m) was mounted on a 
mast above the tower and all lower anemometers were mounted on booms 
into the wind. Studies showed that the boom-mounted anemometers were 
too low by amounts up to 15% because of tower interference. It can 
readily be shown that if the winds and temperatures used in the 
analyses are in error in opposite senses the consequent error in the 
calculated values is significant. This is particularly true in cal- 
culations of the Richardson numbers. 

A study of Figs. 1-4 reveals that the variation of the power 
law exponent p depends on the state of turbulence in the atmosphere. 


In the free convection regime the value of p remains approximately 
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constant while in the forced convection regime p increases rapidly 
with increasing values of Ri. The comparison of Figs. 1-4 shows that 
the discontinuity in the p : Ri curve occurs in the vicinity of 

Ri = -0.03, in agreement with Priestley's (1959) value of the free 
convection-forced convection transition. 

Figs. 5-8 show the comparison of the accuracy with which the 
power law, the logarithmic law, and the log-linear law represent the 
observed data. If the observed wind speeds fitted a particular law 
with one-hundred per cent accuracy, the weighted correlation coef- 
ficient would be equal to unity. It may be observed from the graphs 
that the Kerang data are represented most accurately by the Monin- 
Obukhov log-linear law while the power law is the least valid. For 
the 0.5-16 m layer at Suffield it is observed that the power law 
represents the data as well as or better than either of the other two 
laws with the exception of a short interval in vicinity of Ri = -0.01. 
This finding is in agreement with that found in an independent study 
by Johnson (1959) when he compared the logarithmic and the power law 
using the Suffield data. Figs. 7 and 8 show that in the 2-92 m and 
16-92 m layers the log-linear law represents the data best in the 
near-adiabatic vertical temperature stratification while in the regions 
far from adiabatic the power law gives a better fit. This result is 
in agreement with the theoretical derivation of the Monin-Obukhov 
log-linear law which, according to theory, is valid only in the near- 


adiabatic conditions. 


yvibigs1 esesetoni q 



















sits 


i 
- 


ee ve . rue 
smige1 motjosvnoo beozol oft nb alindw ans3e 


361 ewore A-! .egti to noetisqmoo sit iff Jo seulav gatesstont dalw 
io. yJiclotv ad3 at aavoso avauo FA: q ef3 al ySiunisnoselb adg 
a - 
, 2 r : » es . o r ' 
so7t sls sulav (@¢8l) e'yelseettd Astw inemeetge ni £0.0- = be ; 
; P - Pea ae 
soljianetji noldosvnes beore)-notssevnoo 
= 
* 
in bel f : 3 ab ie eonel 2.2 an td 
ai3 doldw d3iw yastvoos sdi to mosltsqmos 91 wore 8-2 .eglht 
‘ - 
3 Jnsastqst wel wssaii-gol ed3 bos .wal olmitiazerol sit .wel Iswog 7 
> 
‘ei rsivoijisgq & bassit abesge bniw bevaeedo of3 TI .sitsb bevasedo 
=. 
a . r , . ; 4 j . if ' 
-te09 nolisieszo09 beidgiow of2 .yosiv508% tnso 19q besbausi-eno d3iw 
- - 
edigatg of3 mori bevteedo sd yam 3] yjinu o3 Isupe ed bluow jastolt 
; = : ory 
‘ 2 ne 
-ninoM ais yd yiesstssw59% om bsijnesseszgey 9 B3aD gnsyiexX odd jad3 o 
—- 
oF bilsv sess! 3 ai [ xewoq of3 olidw wal zwenil-gol voddudd 
Piers 
a 
wel aswoq edd teni3 hevisado el 31 bislitw2 38 seysi m 31-2.0 afd 
F . ; 
ows tadjo edi to tonite nsd+ tetied ro as Ilaw es stab sid 2309893997 - 
i 
a 
rh Cj "| ve y s r F : wy eT ei 4 r ¢ aA r*eacy afi 3 1: 8 oo 
sibs an 4h oO PILMNLOL rm BVIS3H 230n8 aV NOLIsSG 32329 ati3 Ajiw swe. 7 
eal a 
ybuse 3osbasqebnt os nt bayol jad3 ditw dnemse7gs nt ek gatbat? aidT 
wy. 
wel wewog erid bhi oimdiiusgol sid beiageoo on pesiw (C20!) noandol yi —stié‘éE= 
— C00 ot so Da > ha ; ier Se P Pa - a 
O78 M .k=s BOI Gi 3683 wota G bos \ .egiyi .ageb bleittiu? ed3 gotau 
- “ ia 
7 - a 5s > = £ " » - ~~“. oO a os , . s i 7 
9d3 mi jasd stab arid etmosexqex wel asemif-gol asi: azoyel @ Se-8t 
_ 


enoiget ei3 ait slidw sol 


ai tiuasz eliT 


vodaudO-cinoM eds 


7 


“— sen off2 ak vino biiav ai ,yt0ed3 03 gaibtoozs ,doldw wal teed age 


Poe 
olisdalbs--: 






a 
IB OLLLI67398 STU ISIS¢Mm 


Zz 


Yasied s oljsdaibs moxi 2683 


eovig wel rewog en 









5 wad 
PS q M . 5 7 
to molsavixeb Isotierzoerds eff d3iw Insmes 36 ot 
; : i he ary 
j » 





> | 


a —— 





Ey 

The comparison of accuracy with which the Brooks log-power 
law, the modified log-linear law, and the Monin-Obukhov log-linear 
law represent the data is shown in Figs. 9-12. It may be observed 
that for the Kerang data the ar teereacea between the three laws are 
very small. There exists a systematic tendency for all three laws 
to fit the data more accurately in vicinity of the adiabatic condition. 
However, there is no significant departure of any particular law from 
the remaining two. On the other hand, for the 0.5-16 m layer at 
Suffield the differences among the three laws are negligible for 
Ri <0 while for Ri > 0O the differences become more pronounced. 

It appears that the Brooks log-power law becomes less valid under the 
inversion conditions, whereas the modified log-linear law represents 
the data well in this stability regime. In the Suffield 2-92 m and 
16-92 m layers all three laws appear to fit data equally well. The 
slight differences which do exist lack a systematic pattern and cannot 
be relied upon to determine the superiority of one law over another. 
This result appears to coincide with findings by Bernstein (1966) 

and Hansen (1966) when they analysed several sets of data in terms 

of different forms of the log-linear law. 

The question of how significant are the differences between 
the power law, the logarithmic law, and the Monin-Obukhov log-linear 
law may be answered by studying Figs. 13-16. The vertical line seg- 
ments represent the five per cent level confidence intervals for 
different laws. It may be noted from the diagrams that for Kerang 


data the Monin-Obukhov log-linear law is significantly better than the 


afte 
Tawog-gol exeo78 sis notdw dilw ysatuvoon 


tsenil-go! voriat bns 


bavtaado sd y nwone 


.w8I ssefti-go0l beltibom ons . 


: i 


to soelveqmos odT - 





















a 


ss 


+ fal i 
<é 


ai eiab art tnseetqs7 § 
| 7. 


- 
5% sd? —= ‘Jat 








" r 
978 ews! sexni sit asswied esonsysitib add sisb gnats + 
ewsfl sexda Ifs 1ot yonebnsd oltsmadeye s agalxe stefiT .ilame visv AS 
7” 
a , 
M0isibnos sitedsibs efi to yiinifotv at ylsdsatvcon oom |3eb ed 332 os 
7 
mom? wei rsfuoltiag wre to sivtyaqeb saBoti ins on @f etsdi3 .tavewoH 
38 tsysi m of-c.0 aft rol , bred deiito sd3 nO .ow2 gniniemet ods 
a 
1602 sid} on otf ewel ssid ot gnome esome7st?rb aris biertive — 
» 3 YY 7 i i ia : 
s9nvOnoIG sx smoos¢d easomstsitib sri 0 tf +wet slidw O> 19 
a. 
7 
3 x9bnu Stlav easl esmoossd wel xswhg-golf exioors sit sand axBsqqs 3I 
~ Pa 
2irsestgs1 wal + i-gol Beiitibom sis asetedw .s Jibsos mofersvar 
, - 
a c e: £ » eT on . > re. | ta ey - ’ + 
bas’ m $@- £ u2 eft al migor yiilideje efda ot ilew sisb sd 
wit [low vifsups sts! OJ TSsggs ews! soxis lis ete m S@-87 
- Sey ie “ 4 mse ° 4 . 4 fon bd ; 
Jonnss bos nrsdgsq ottsmeteve aosi jetxs ob dotdw eaons7s32ib Srigtle 
x 
<tetdons tsyo waf smo to yvilzols é 33 snimisisb o3 noqu belle® bed - 
i a > ~ - oe I . J <a cg 64 
(OogeJ risajenrwds yd rt Brtt ditw abionios of aIssq tiuasy at - 
ot 2 _ ~ « fs on ol r = on al ot a fr os 7 
ents2 Of 6385 to ejsa larevee boeyvlens yads nedw (0901) nsana thie 
.w8l rwenti-g0[ sfd3 to amrol tnexs?? 
_ 
MeswIsd ason 


wesli-gol voravudO-minoM sd1 bas ,wel sims 
-g9a ontl [aot3asy sit 


402 elevietat sonsbtisos eval 3n99 199g 


ay x02 









ot wlanenttinde 


a 





wat Ban 
PAS 


1p 


e79711b sd2 918 jnsviitngie wod to molseeup sdiT 


Gi-EL .egtt enivbise te betawans ed ea ar 


jed3 emezgsib sy mor besten od ee aa! 





Hae 







= >) 
irtsg01 ofa wel on ae 












‘evi orld snovexgor snl 


oi 
rn 


odslede 
7 ay 





ayan 

power law in vicinity of Ri = -0.03 while in the other regions the 
confidence intervals overlap. For the Suffield data the significant 
differences occur in the 2-92 m layer between the Monin-Obukhov log- 
linear law and the logarithmic law in the interval 0.00< Ri< 0.20 
and between the power law and the logarithmic law in vicinity of 

Ri = -0.20 and Ri = 0.40. 

Figs. 17-20 show the confidence intervals for the Brooks log- 
power law, the modified log-linear law, and the Monin-Obukhov log- 
linear law. It appears that the confidence intervals for these three 
laws overlap for nearly all stability conditions, This result suggests 
that the differences among the two forms of the log-linear law and the 
log-power law are not significant at the five per cent level. 

The relatively wide confidence intervals appearing in Figs. 
17-20 may be somewhat misleading because, as shown in Appendix B, the 
width of the confidence intervals depends on the numbers of observations. 
This fact would suggest that the confidence intervals calculated not 
for each subgroup, but rather for each set of data may provide a better 
estimate of the significant differences between the different laws. 
Figs. 21-24 show the average weighted correlation coefficient, together 
with the corresponding confidence intervals, for each set of data. 

It may be observed from Fig. 21 that for the Kerang data the log-linear 
law is significantly better than either the logarithmic law or the 
power law. Whether the logarithmic law is significantly better than 
the power law is open to question because the five per cent level 


confidence intervals for these two laws just touch one another. 
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For the Suffield 0.5-16 m layer the power law is significantly 
better than the logarithmic law. However, there is no significant 
difference between these two laws and the log-linear law. In the 
Suffield 2-92 m layer the log-linear law is significantly better than 
the logarithmic law while there is no significant difference between 
the power law and either of the other two laws. Lastly, in the Suffield 
16-92 m layer confidence limits from all three laws overlap, hence 
no one of the laws is significantly different from the other two. 

The comparison among the average weighted correlation coef- 
ficients and the corresponding confidence intervals, for Brooks log- 
power law, modified log-linear law, and Monin-Obukhov log-linear law 
is indicated in Table 5.1. A study of Table 5.1 reveals that there 
are no significant differences among these three laws, a confirmation 


of the results shown in Figs. 17-20. 
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CHAPTER VII 
CONCLUS IONS 


The study of the wind profile measurements obtained at 
Suffield, Alberta, and at Kerang, Australia, shows that the value 
of the power law exponent p is smallest during slight and moderate 
superadiabatic conditions but increases slightly in extreme lapse 
conditions. Moreover, relatively large values of p exist during 
inversions. The numerical values of p obtained in this study are 
in agreement with those obtained at other locations by different in- 
vestigators. It was found that the value of p _ remains approximately 
constant in the free convection regime and it increases rapidly in the 
forced convection regime. 

The values of the parameters u, and z4 appearing in the 
logarithmic and log-linear law, were evaluated from the Suffield wind 
profile data using the graphical method. The average values of 


u and Zo for Suffield terrain were estimated to be 0.34 m Pose 


* 
and 3 cm respectively. In addition, the parameters uU,, 25 and 

D were calculated using Lettau's method. According to this method 

the average values of u,, z, and D, for Suffield terrain, are 

0.34 m seat: 3 1.5 cm and -3.2 cm respectively. Because of the nature 
of the vegetation cover at Suffield the value of Zo given by Lettau's 
method is considered more accurate than that given by the graphical 
method. 


The parameters a, Y' and Ly), which appear in the Monin- 


Obukhov log-linear law, the Brooks log-power law, and the modified 
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log-linear law respectively, were calculated from the wind profile 
data. It was found that the value of the Monin-Obukhov constant o 
varies with Richardson number and also with the height of the layer 
considered. On the other hand, the parameter Y' was found to be 
independent of the height range but it varies with location and with 
stability. The length L, was found to assume large values under 
unstable conditions and small values under inversions. 

The comparisons among the power law, the logarithmic law, and 
the Monin-Obukhov log-linear law reveal that in the near-adiabatic 
conditions the log-linear law represents the observed data more 
accurately than either of the other two laws. However, when the 
conditions are not near-adiabatic the power law gives a better rep- 
resentation, although the difference is not significant at the five 
per cent level. The logarithmic law which is a special case of the 
log-linear law for Ri =0, is found to be valid only under adiabatic 
conditions. 

The comparisons among the Monin-Obukhov log-linear law, the 
Brooks log-power law, and the modified log-linear law show that the 
differences among these three laws are very slight and are not 
significant at the five per cent level. It would seem that the question, 
which of these three laws represents the data best, cannot be settled 
from the data which were available for this study. 

Because of the tower interference in the Suffield 0.5-16 m 
layer wind measurements, discussed in Chapter VI, the results obtained 


for that layer are considered to be less significant than those 
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obtained for the other layers. 

It may be observed from the graphs that the Kerang 
measurements display less scatter than do the Suffield data. This 
fact would suggest that a high degree of accuracy is required in the 
data to distinguish among the different wind profile models. There- 
fore, in the future experiments considerable care should be taken to 


ensure that a high degree of accuracy is achieved. 
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APPENDIX A 


SUFFIELD WIND PROFILES FROM WHICH u, AND Zo 


WERE OBTAINED 
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TRIAL 2u 

DATE 1 /2 / 60 

TIME 1510 MST. 

SAMPLING PERIOO 60 MIN. 

WNO,OIR. 170 OEG. 

TG (M-1/2M) =-0.3 OEG, (F) 
TG (13M-1M) =MSG. DEG. (F) 
TG (92M-1M) = MSG. DEG. (F) 
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TRIAL 25 

DATE 8 2 / 60 

TIME 1641 MST. 

SAMPLING PERICO 36 MIN. 

WNO.OIR. 70 OEG. 

TG (UM-1/2M) =0.0 OEG. (FI 
TG (13M-1M)=MSG. OEG. (F) 
TG (92M-1M) = MSG. DEG. (F) 
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TRIAL 26 

DATE i2/ / 60 

TIME 1156 MST. 

SAMPLING PERIOO 60 MIN. 

WND.OIR. 195 OEG. 

TG (UM-1/2M) = -0.4 DEG. (F) 
TG (13M-1M) = MSG. OEG. (F) 
TG (S2M-1M) =MSG. OEG. (F) 
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TRIAL 32 

DATE 23/7 / &3 

TIME 1857 MST. 

SAMPLING PERIOD 60 MIN. 

WNO.OIR. 265 O&G. 

TG (UM-1/2M) =-0.3 DEG. (F) 
TG (13M-1M)=MSG. DEG. (F) 
TG (92M-1M) =MSG. DEG. (F) 
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16.00 20.00 


24.00 






Var 04 3I1A0A9 QUIN 
Ss .314 





HEIGHT (METERS) 


0.00 


4.00 


8.00 
U 


—90-— 


WIND PROFILE NO. 
fiGs 90 


12.00 
(MPH. } 


Excl 


TRIAL 36 
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TIME 1526 MST. 
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TG (UM-1/2M) =-0.3 DEG. (F) 
TG (19M-1M) =MSG. OEG. (F) 
TG (92M-1M) = MSG. OEG. (F) 
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APPENDIX B 
STATISTICAL FORMULAE 


I. Simple linear regression 
Given n sets of observations (z;, Uj), i=1, 2, ... Qn, 
of an independent variable z and a dependent variable U_ then a 


straight line in the form: 
U = “az + Db (Be) 


can be fitted to these data by the method of least squares. The sample 


mean and the standard deviation of z are given by: 


zhes = nied. (B.2) 





and 
ah ge B(z, - 2)° (B. 3) 


Zz n- l 


respectively. Similar expressions hold for U. The least-squares 


estimates of a and b in Eq. (B.1) are given as: 


Ta(2 mec ae lua U) 


ee eee (B.4) 


and i 
b = U- az (B.5) 


The sample correlation coefficient between z and U is: 


as, 
(B.6) 
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2s Fisher's z' - statistic 

According to Brooks and Carruthers (1953), if the values of 
a correlation coefficient (r), computed from two or more rather short 
sets of data, are known a better value of r may be obtained by 
combining them. This is true only if it is reasonable to suppose that 
the different sets of data are samples from the same population. The 
correlation coefficients are combined by transforming them first into 


Fisher's z' - statistic given by: 


gi = - [1n (ee) sine 1-9 1) | (Bat) 


This statistic has a nearly normal distribution and a standard error 


equal to: 





Oo 1 
é (B.8) 
(n - 3)% 


Next, a weighted mean of the z ’s is found, the weights being 


proportional to n - 3 The weighted mean of z'’s is converted back 


to xr by equation: 


exp Cay en" 
explie2 |) Peel Sate. 


To calculate the 5% level confidence limits the values of 
z' +1.960 and yl = 1 59610-. “ane calculated first and then transformed 
to xr by Eq. (B.9). It may be noted that the relative size of con- 


fidence limits is dependent on the number of observations. 
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APPENDIX C 


SELECTED WIND PROFILES FOR 


SUFFIELD 2 - 92 m LAYER 
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TRIAL 1 

DATE 7 /10/ 63 

TIME 1603 MST. 

SAMPLING PERIOD 10 MIN. 

WNO.OIR. 275 DEG. 

TG (UM-1/2M) =-1.4 DEG. (F) 
TG (13M-1M) =-1.0 DEG. (F) 
TG (92M-1M)=2.1 DEG. (F) 
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TRIAL 1 

DATE 7 AO/ 63 

TIME 1613 MST. 

SAMPLING PERIOD 10 MIN. 

WND.DIR. 265 DEG. 

TG (UM-1/2M) =-0.8 DEG. (F) 
TG (13M-1M) =-1.0 DEG. (F) 
TG (92M-1M)=0.2 DEG. (F) 
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DATE 7 AO/ 63 

TIME 1623 MST. 

SAMPLING PERIOO 10 MIN. 

WND.OIR. 276 DEG. 

TG (UM-1/2M) =-0.8 DEG. (F) 
TG (13M-1M) =-0.4 DEG. (F) 
TG (92M-1M) =-0.1 DEG. (F) 
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